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Abstract 

We show that C-odd gluelumps can be successfully described as bound states of a single trans- 
verse constituent gluon evolving in the flux-tube-like potential generated by a static color-octet 
source. The use of a helicity degree of freedom rather than a spin one for the constituent gluon for- 
bids the states that are not observed in lattice QCD. Our model leads to a gluelump mass spectrum 
in remarkable agreement with the available lattice data provided that an additional parity-splitting 
mass term is introduced. We argue that such a term is due to instanton-induced interactions in 
gluelumps. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) allows the existence of pure gauge states as glueballs 
and gluelumps. Glueballs are bound states of the gluonic field alone, while gluelumps, to 
which this work is devoted, are bound states of the gluonic field plus a static color-octet 
source. A rdative gluelump mass spectrum has beeu first computed in lattice QCD Q, 
while the determination of the absolute masses has been achieved in ref . |2J . Gluelumps were 
originally a first attempt to model gluon-gluino bound states, allowed by supersymmetric 
theories But it has been soon realized that gluelumps are of importance in hadronic 
physics too. Indeed, if the static source is seen as a pointlike heavy quark- ant iquark (qq) 
pair, the gluelump mass can be interpreted as the energy of a static qq pair in the limit 
where their separation (r gf ) vanishes. Such an energy has been computed in lattice QCD 
for r q q ^ [sj], and it indeed appears that the extrapolations of the different energy levels at 



r qq 



match very well with the gluelump masses [2j. The study of gluelumps is therefore 



in which there will be an 



strongly linked to the understanding of heavy hybrid mesons [3|, 
increasing experimental interest with future experiments like BESIII, GLUEX and PANDA. 

Apart from lattice QCD, gluelumps have been studied within different approaches: Bag 
model [5I, Coulomb gauge QCD |6|, pNRQCD [4], and potential models [7]. In the framework 
of potential models, gluelumps are usually seen as a single constituent gluon evolving in 
the potential generated by the static source. The constituent gluon is generally assumed 
to have a spin degree of freedom. However, it has been recently shown that the lattice 
QCD data concerning pure gauge hadrons were rather compatible at a qualitative level with 
bound states of massless, helicity-1 (transverse), constituent gluons |8|. Moreover, several 
arguments are given in this last reference indicating that a potential model is a relevant 
approach to study gluonic hadrons. By helicity, we mean that the zero projection of the 
gluon spin is forbidden. It has already been shown that describing the C-even glueballs as 
bound states of two massless helicity-1 gluons interacting via a flux-tube-inspired potential 
leads to a good agreement with lattice QCD In particular, the use of helicity as intrinsic 
degree of freedom imposes constraints on the allowed glueball quantum numbers: Light 
vector glueballs, that are not seen in lattice QCD, are forbidden. This is not the case with 
spin-1 gluons. 

Motivated by the above discussion, we propose in this paper to reconsider the gluelump 
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spectrum. We first discuss the gluelump quantum states, referred to as gluelump helicity 
states, in sec. [Til We then propose a semirelativistic Hamiltonian, based on the flux tube 
model, in sec. IIHI The C-odd gluelump spectrum is numerically computed and discussed in 
sec. IIVI while a summary of our results is given in sec. |Vj 



II. HELICITY STATES FOR GLUELUMPS 



The lowest- lying gluelumps have a negative charge conjugation and are lighter than the 

lowest-lying glueballs. We have shown in ref. 8j that this fact was compatible with an 

interpretation of gluelumps as bound states of a single constituent gluon within the potential 

generated by the static source, while glueballs should be made of at least two constituent 

gluons. Notice that the static color-octet source has the vacuum quantum numbers ++ . 

Our basic assumption is that a constituent gluon is a massless transverse particle. The 

quantum states describing bound states of such particles can be formulated within the 
□ 

helicity formalism [10[. We refer the reader to this last reference for a detailed description 
of this formalism. It is enough for our purpose to recall that, within the helicity formalism, 
there are two families of one-gluon states. Explicitly written in a standard | 25+1 Lj) basis, 
they read [8j 



1 |3 J+U 



2J + 1 

with P = (-) J , (la) 
\B;(J>l) p -) = -\ 3 Jj) 

with P = (-) J+1 . (lb) 

The use of helicity imposes that J > 1; such a constraint is not present with a spin-1 gluons. 
The fact that J = states are forbidden is in agreement with lattice QCD since no light P ~ 
gluelump is observed [l|. Moreover, the lightest gluelumps are l p ~ ones as expected from 
eqs. (CQ). Notice that the helicity formalism also reduces the number of possible gluelump 
states, and thus the arbitrary of the model. For example, \ A\ 1 ) is the unique 1 gluelump 
helicity state, while both | 3 Si) and | 3 -Di) would be allowed with spin-1 gluons. 

Thanks to the above decompositions, simple computations show that the matrix elements 
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of the different operators usually appearing in potential models (square orbital angular 
momentum, spin-orbit, tensor, etc.) have the same values for both families. One has indeed 
(L 2 ) = J{J + 1), (S 2 ) = 2, (L ■ S) = -1, and (S 2 - 3(5 • rf) = -1, where the average 
values are computed with either the state ( {Taj) or the state ( llbl) . 

Gluelumps with a positive charge conjugation cannot be one-gluon states: At least two 
gluons are needed jf], [sj] . This is coherent with the fact that the currently known C-even 
gluelumps (0 ++ and 1 h ) are well heavier than the C-odd ones. The corresponding helicity 
states can be obtained by coupling the two constituent gluons to a given total spin j and 
parity tt, and by recoupling these two gluons, seen as a massive particle of spin j and parity 
7r, to the static source following a general procedure described in ref. Ill ] . General two-gluon 
states are extensively studied in ref. [9]. Here it is sufficient to recall that the lightest totally 
symmetric two-gluon configuration is reached for j n = + . The two-gluon cluster and the 
static source, being both massive objects, can be directly coupled in the usual | 2S+1 Lj) 
basis. The two lightest states should be mainly the | Sq) and | Pi) ones, which have the 
quantum numbers ++ and 1 h respectively, in agreement with lattice QCD. Notice that in 
the following we focus on C-odd gluelumps only. 



III. THE MODEL 



At the dominant order, the flux tube Hamiltonian describing a gluelump made of one 
massless constituent gluon reads 

H = + , (2) 

r 

where p is the momentum of the constituent gluon and where r is the separation between 
the constituent gluon and the static source. Notice that r = \r |. This Hamiltonian is spin- 
independent but the gluon helicity is already included at the level of the gluelump helicity 
states (00). The formalism is not explicitly covariant, but the semirelativistic kinetic term 
^/fi 2 is relevant to describe the dynamics of a massless particle, a is the energy density of the 
straight flux tube linking the constituent gluon to the static source. Assuming the Casimir 
scaling hypothesis, one can write a = (9/4) a, with a the energy density of a fundamental 
flux tube, that is the flux tube in a meson. The flux tube stands for the nonperturbative 
part of the QCD interactions and generates the confinement. It is a dynamical object that 
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carries both energy and angular momentum but, at the lowest order, it reduces to a linearly 
rising confining potential. The flux tube clearly dominates the long range interactions. 
Nevertheless, other mechanisms come into play at short distances, i.e. one gluon exchange 
processes. The Coulomb part, where a s is the strong coupling constant, is the dominant par t 
of the one gluon exchange potential between the constituent gluon and the static source 12]. 

Lattice QCD mass spectra are generally expressed in units of a fundamental length scale r 
that can be linked to a by the relation a = 1/r 2 .. Most of the final lattice error bars actually 
comes from the determination of this length scale r , and it is thus worth comparing our 
model with the dimensionless lattice QCD spectrum r^M^. To this aim, it is convenient 
to work with the new coordinates q = r$p and x = v/tq. Then, Hamiltonian (TSJ) becomes 

9 a s 

r H = a/? 2 + -x - 3— . (3) 

4 x 

This rescaling eliminates a free parameter in our model: Only a s has now to be fitted on 
the data. Before making explicit numerical computations, it should be written that 

This particular value of the square orbital angular momentum is a consequence of the helicity 
states ([T]), as we already mentioned in the previous section. 

The Hamiltonian only depends on J through eq. (j4l . It means that the gluelumps 
of quantum numbers J ±_ are degenerate if only H is used to compute the gluelump mass 
spectrum, in obvious contradiction with lattice QCD. Adding the usual spin- dependent 
corrections of Fermi-Breit type would not lift this degeneracy since, as shown in sec. HIl 
the matrix elements of the spin-spin, spin-orbit, and tensor operators are identical for both 
gluelump helicity states. This leads us to introduce a phenomenological mass term of the 
form 

r AM = (-) J r o7 P, (5) 

splitting the gluelumps with opposite parity but with the same total spin J. A priori, 7 could 
be a complicated function of J, x, etc. Nevertheless, we assume here that 7 is constant and 
positive. The meaning of this extra term will be discussed in the next section. Notice that 
the need for an interaction term going beyond the standard pair interactions in gluelumps 
has also been pointed out in the Coulomb gauge study of ref. [6(. 
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TABLE I: Gluelump spectrum computed from eqs. 
ref. J 



compared with the lattice QCD data of 





r M Lat [2] 


r M 


1+- 


2.25(39) 


2.25 


1~ 


3.18(41) 


3.19 


2~ 


3.69(42) 


3.70 


2+- 


4.72(48) 


4.64 


3+- 


4.72(45) 


4.70 


4~ 


5.41(46) 


5.52 


3~ 




5.64 



IV. RESULTS AND DISCUSSION 



The mass spectrum of Hamiltonian 
using the Lagrange n 
the total gluelump mass 



, denoted as r M , can be accurately computed 
by using the Lagrange mesh method [L3J. Adding the mass term (j5]), one can finally obtain 



r M = r M + r AM, 



(6a) 



where r$ Mq is the eigenvalue of Hamiltonian ([3]). 

As it can be seen in table [U a remarkable agreement between our model and the available 
lattice QCD data is reached provided that the following values are used for the parameters 



O's 



0.38, r 7 = 0.47. 



(6b) 



Note that a s ~ 0.4 is a standard value in light hadron physics. In particular, the light 
meson [yj] and glueball 9] mass spectra can be satisfactory reproduced within the flux tube 
model by using such a value for the strong coupling constant. Notice that the results of 
tabled] can be converted in physical units with the typical conversion factor Tq 1 = 440 MeV. 
One has then 7 m 0.2 GeV. 

Without the phenomenological term (jHJ), the masses of gluelumps with the same total 
spin but opposite parity would be equal, in disagreement with lattice QCD. Once it is 
introduced, it first reproduces the correct ordering of the 1 ±_ and 2 T_ states. Second, it 
lowers the 3 +_ mass so that it is nearly degenerate with the 2 H , as observed in ref. 
Finally, it causes the 4 gluelump to be lighter than the 3 one. This is not an intuitive 



TABLE II: Mass gaps (171) computed from eqs. 

n 

ref. [l[. 



and compared with the lattice QCD data of 





roM* at [1] 


r M R 


1— 


0.933(18) 


0.94 


2— 


1.438(25) 


1.45 


2+- 


2.467(92) 


2.39 


3+- 


2.468(60) 


2.45 



result since Hamiltonian ([3]) alone states that gluelumps with higher J should be always 
heavier than states with lower J. 



The results of our model can a 
that has been computed in ref. 



so be compared with the relative gluelump mass spectrum 
I]. In this last reference, the mass of the l +_ gluelump 
(Mi+-) is not known, and only the mass gaps between the various gluelump states and the 
l + ~ one can be computed. But, the accuracy of the results is better than in ref. |2|. We 
denote 

r M R = r M - r M 1+ - (7) 

the mass gaps. As it can be observed in table [Til our predicted values compare very favorably 
with lattice QCD. In particular, the mass gaps concerning the 2 and 3 H gluelumps are 
independent of the term (J5j). The fact that they are compatible with the lattice results is a 
good check of the relevance of Hamiltonian ([3]) alone. 

Although leading to nice results, the mass term (jSJ) is introduced ad hoc. Obviously, its 
origin cannot be linked to some short-range correction of one gluon exchange type since all 
the operators appearing in such interactions only depend on J, not on P. One could think 
about higher-order corrections, with operators such as L A for example, but these corrections 
are not proportional to P in general. Consequently, this term should be nonperturbative, 
and instanton-induced forces are a relevant candidate. Indeed, as a consequence of the self- 
duality of the instanton's field strength, it has been argued that the instanton contribution 
is generally equal in magnitude but opposite in sign for hadrons with the same total spin and 
charge conjugation, but with opposite parity [15], just as the mass term (jSJ) is. Although 
instanton-induced interactions are well known for mesons, it is not the case for gluonic 
hadrons. From refs. [jij, it can be deduced that instantons induce a mass term of the form 
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—IPSjfi for C-even glueballs. Such a term has been successfully used in ref. which is 
clearly a two-gluon generalization of the present work. We are thus led to conjecture that 
our phenomenological term (jSJ) originates in instanton-induced forces in gluelumps. To our 
knowledge, no study of these interactions has been performed yet and this is a task out of 
the scope of this paper. We hope that such results will be available in the future in order 
to check if the form (—) Jr yP eventually emerges as an effective instanton mass term in 
gluelumps. 

The (— ) J factor in eq. (0) could seem rather surprising. It is actually such that, for a 
given J, the helicity state ( llbl) is always lighter than the state ( ITal) . The \B) family, the 
lightest one, corresponds to the gluelump states of natural parity^that is r\ g (— 1) J , rj g being 
the intrinsic parity of a gluon. We already pointed out in ref. |8J| that the gluonic states 
appearing in lattice QCD (glueballs and gluelumps) are always such that those with natural 
parity are lighter. Our mass term (j^J), as well as the currently known gluelump spectrum, 
is in agreement with this fact. 



V. CONCLUSION 



In summary, we have described the C-odd gluelumps as bound states of one transverse 
gluon within the potential generated by a static color-octet source. The gluelump quantum 
states are computed thanks to the helicity formalism, and the Hamiltonian corresponds to 
the dominant order of the flux tube model with one-gluon exchange potential. A very good 
agreement with lattice QCD is reached provided that an extra term (—) J/ yP is introduced. 
Although further research in that domain is needed, we expect that this term is intrinsically 
nonperturbative and comes from instanton-induced interactions in gluelumps. 



Acknowledgments 

The author thanks Claude Semay for valuable discussions about the present work. 



[1] M. Foster and C. Michael, Phys. Rev. D 59, 094509 (1999) [hep-lat/98li010] . 
[2] G. S. Bali and A. Pineda, Phys. Rev. D 69, 094001 (2004) [hep-ph/03101~30j . 



S 



[3] K. J. Juge, J. Kuti and C. J. Morningstar, Nucl. Phys. Proc. Suppl. 63, 326 (1998) 
[hep-lat/970 9131|; Phys. Rev. Lett. 90, 161601 (2003) [hep-lat/0207004] . 

[4] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B 566, 275 (2000) 
|hep-ph/9907240| . 

[5] G. Karl and J. E. Paton, Phys. Rev. D 60, 034015 (1999) l; hep-ph/9904407] . 
[6] P. Guo et al, Phys. Rev. D 77, 056005 (2008) [arXiv:0707.3l"56| . 



[7] Yu. A. Simonov, Nucl. Phys. B 592, 350 (2001) |hep-ph/0003114|; E. Abreu and P. Bicudo, J. 



Phys. G 34, 195207 (2007) [hep-ph/0508281 1; V. Mathieu, C. Semay and F. Brau, Eur. Phys. 



J. A 27, 225 (2006) [hep-ph/0511210] . 

[8] N. Boulanger, F. Buisseret, V. Mathieu and C. Semay, EIII\E0806. 3174. 

[9] V. Mathieu, F. Buisseret and C. Semay, Phys. Rev. D 77, 114022 (2008) [arXiv:0802.0088] 
[10] M. Jacob and G. C. Wick, Ann. Phys. 7, 404 (1959). 
[11] G. C. Wick, Ann. Phys. 18, 65 (1962). 

[12] V. Mathieu and F. Buisseret, J. Phys. G 35, 025006 (2008) [hep-ph/0702226| . 

[13] C. Semay, D. Baye, M. Hesse and B. Silvestre-Brac, Phys. Rev. E 64, 016703 (2001). 

[14] F. Buisseret, Phys. Rev. C 76, 025206 (2007) [arXiv: 0705 .0916] . 



[15] N. Kochelev and D. P. Min, Phys. Lett. B 633, 283 (2006) [hep-ph/ 0508288|. 
[16] T. Schafer and E. V. Shuryak, Phys. Rev. Lett. 75, 1707 (1995) (hep-ph/9410372| ; H. Forkel, 
Phys. Rev. D 71, 054008 (2005) |hep-ph/0312049] . 



9 



